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The vortex glass state formed by magnetic flux lines in a type-II superconductor is shown to 
possess non-trivial three-body correlations. While such correlations are usually difficult to measure 
in glassy systems, the magnetic fields associated with the ffux vortices allow us to probe these via 
muon-spin rotation measurements of the local field distribution. We show via numerical simulations 
and analytic calculations that these observations provide detailed microscopic insight into the local 
order of the vortex glass and more generally validate a theoretical framework for correlations in 
glassy systems. 



PACS numbers: 74.25.Qt, 76.75.-l-i, 61.20.Gy 

In systems which possess long range order, such as 
atomic crystals, the local arrangement of particles is eas- 
ily obtained from scattering experiments. In disordered 
systems, the average correlation between the positions of 
two particles can be measured by scattering techniques, 
but inferring anything more about the local geometry 
is a far more subtle issue. Little is known experimen- 
tally about correlation functions of higher order. Mea- 
surements of three-body correlation functions for colloids 
imaged in a quasi-two-dimensional geometry have been 
reported recently [H. However, hulk measurements of 
three-body correlation functions in any system are still 
unavailable and our understanding of such correlations 
derives mainly from simulations. This Letter reports a 
study of the local structure of the vortex glass phase in 
a bulk type-II superconductor. The vortex glass phase 
provides an example of a glassy system where the local 
geometry in the hulk is uniquely amenable to investiga- 
tion, due to the magnetic fields associated with the vor- 
tices, which we measure by the muon-spin rotation (/zSR) 
technique. By coupling these measurements with Monte 
Carlo simulations and analytic calculations, we demon- 
strate both the existence of non-trivial three-body corre- 
lations in the flux-line array and a theoretical framework 
in which they may be understood. 

In the mixed state of a type-II superconductor, an ap- 
plied magnetic field penetrates as lines of magnetic flux, 
quantized in units of the flux quantum $0 — h/2e. Such 
vortex lines would form an Abrikosov flux lattice at low 
temperatures in the absence of quenched disorder. As 
the temperature or the strength of disorder is increased, 
ordered arrangements of vortex lines yield to disordered 
ones [1, Q . Weak quenched disorder converts the crystal 
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FIG. 1: (color online) The field distribution function n{B) 
for two values of the applied magnetic field H at T=5 K: (a) 
H = 40 mT and (b) H = 0.5 T. The insets show the same 
data, plotted over the same range, with the y-axis plotted on 
a logarithmic scale, indicating the errors on the points. The 
curves are normalized at the mode of the distribution, Bpi,. 



into a "Bragg glass" with quasi-long range order in trans- 
lational correlations [1]. At stronger disorder, "vortex 
glass" states with short-ranged correlations are obtained. 
Neutron scattering measurements support the proposal 
of a power-law decay of translational correlations in the 
Bragg glass phase [5j. In contrast, structure and correla- 
tions in vortex glasses remain little understood. 

Our experimental system is Lai.9Sr.iCu04_5 (LSCO), 
an underdoped high-Tc superconductor with properties 
which amplify the effects of thermal fluctuations and 
quenched disorder. It was recently shown using /iSR 
measurements, on the same high quality crystal as used 
in this experiment, that there is a field induced transition 
in LSCO to a vortex glass phase @. Our present /xSR 
experiments were performed using the GPS spectrom- 
eter at PSI, Switzerland with the field nearly parallel 
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FIG. 2: (color online) The lineshape anisotropy ratio a as ob- 
tained in (i) experiments at 5K (black circles), (ii) calculated 
from Monte Carlo simulations of disordered structures (red 
crosses) and from (iii) a theoretical description in terms of 
liquid state theory (blue squares connected by a solid line) . 



to the c-axis of the crystal. The experimental arrange- 
ment was as described in ref. 0. The novelty of the 
present work is that by relating the third moment of the 
magnetic field distribution to an integral over a three- 
particle structure factor, we are able to provide informa- 
tion about three-body correlations in both the ordered 
and the glassy phases. 

In a /iSR experiment, the probability distribution n(B) 
of the spatially varying magnetic field is inferred from 
the muon precession signal (3]. This distribution re- 
flects the arrangement of vortex lines. We measure n{B) 
as a function of external magnetic field and tempera- 
ture, calculating the second moment ([AS]^), the third 
moment {[AB]^), and the related dimensionless line- 
shape anisotropy ratio a = {[AB]^y/^/{[AB]'^)^/'^, with 
the k*^ moment defined by {[AB]'') = Y.n{B,){B^ - 
{B)f /Y,n(B,). Figs. 1(a) and (b) show field distri- 
butions for two values of the applied magnetic field. As 
discussed in Ref. Q, the data at 40 mT (Fig. 1(a)) show 
a tail on the high-field side and hence positive sign of a 
expected for a lattice (or Bragg glass) structure. How- 
ever, the 0.5 T data (Fig. 1(b)) show a broader and more 
symmetric distribution in which this tail is absent and 
the lineshape at this field is slightly skewed the opposite 
way, signalling that the flux-line structure is not ordered 
like a lattice. Instead, it is in a vortex glass state, which 
dominates the phase diagram in the inset to Fig. 3. 

Fig. 2 shows the behavior of the lineshape anisotropy 
ratio a as a function of magnetic field at 5 K after field 
cooling; all points lie within the vortex glass phase @ 
(inset Fig. 3). The experimental data are the sequence 
of black circles; these data are modeled theoretically (see 
below) by the sequence of red crosses and blue connected 
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FIG. 3: (color online) Regimes in H-T space associated with 
a fixed sign of a = ([AB]3)l/V([AB]^)^/^ and thus of the 
third moment. The inset shows the schematic phase diagram 
for this system in H-T space, obtained from a combination 
of muon-spin rotation, magnetic and Small Angle Neutron 
Scattering measurements ^ . The change of sign of a reported 
here occurs deep within the vortex glass phase shown in the 
inset (dotted line in inset) and over a fairly short range in 
field and temperature. 



squares. Note the reduction in a beginning at rela- 
tively low field values, the precipitous change of sign at 
B ^ 0.35 T, followed by saturation at an approximately 
constant negative value. Fig. 3 shows the variation of a 
over the H ~T plane, further illustrating a change of sign 
from positive to negative values that occurs deep within 
the vortex glass phase. We note that a negative third 
moment of the field distribution is also observed in the 
vortex liquid regime; this is an outstanding problem for 
theories of vortex line correlations jH, Q . 

For a system of flux-lines taken to be rigid along the 
z-axis (field direction) the third moment is obtained from 

{[ABf)^ jj dqidq25(3)(qi,q2)6(qi)6(q2)6(-qi-q2). 



where S'(3)(qi,q2) = ^ (V(qi)^p(q2)<5p(-qi - q2)) [i 
is the triplet structure factor and the proportionality con- 
stant is i?/(27r)''$o- Here, &(q) is the field of a sin- 
gle vortex in Fourier space, while <5/c(q) is the Fourier 
component at wave-vector q of the deviation of the 
flux-line density from its average value. The second 
moment is related similarly to the two-particle correla- 
tion function, i.e. to the conventional structure factor 
S{q) = ji{^p{<l)Sp{—q)). In these expressions, N is the 
total number of vortices and the brackets (• ■ • ) denote an 
average over the sample [ll, 12|. Note that S{q) > 0. 
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The simple London model gives 6(q) = 5/(1 + A^g^^), 
where A is the a6— plane penetration depth and q = 
(q_L)9z = 0). However, this yields a magnetic induc- 
tion which diverges at the vortex core. This unphysical 
divergence is eliminated by multiplying by a "form fac- 
tor" f{q) All the available analytic expressions for 
f{q) give positive values for all q For a perfect tri- 
angular lattice, Eq. 1 reduces to a sum over two sets of 
reciprocal lattice vectors. Each term in such a sum is 
manifestly positive. The form factor reduces the value of 
the third moment, but it cannot change its sign. 

We now relax the requirement of a lattice struc- 
ture. One possibility is to assume uncorrelated lines, i.e. 
S{q) = 1, S'(3Hqi,q2) = 1 This limit yields anoma- 
lously large field (Bq) dependent values for the second 
moment (([AB]^) = Bo^q/AttX^ , in comparison to the 
perfect lattice, for which {[AB]^) = 3.71 x IQ-^^l/X'^). 
It also leads to a positive third moment. Alterna- 
tively, if one applies the "convolution approximation" 
from the fluid literature [lOj . one can express 

5(3) 

m 

terms of the two-particle correlator only: S'^^^(qi,q2) — 
S'(qi)5((72)'5'(|qi — q2|)- However, this approach also 
gives positive {[AB]^) because the integrand in Eqn. 1 is 
manifestly positive. The effects of line wandering can be 
incorporated using expressions for S{q±, qz) obtained via 
the boson approximation Q - numerical values of the mo- 
ments are reduced but again the sign of {[AB]^) should 
remain positive [l^ . 

Negative third moments in vortex glass phases which 
lack long range order can thus only result from vortex 
line arrangements with non-trivial three-particle correla- 
tions beyond the convolution approximation. To validate 
these ideas, we have generated simulations of disordered 
states with these attributes. Simulated annealing tech- 
niques were used to find the ground state of a system 
of 6400 particles (rigid vortex lines) interacting through 
a potential of the form jKQ{r/X), and with a quenched 
one-body potential arising from 3700 randomly placed 
pinning sites. The pinning sites were modeled as at- 
tractive Gaussian wells of the form Vd{r) — — j'e~^ ^^p, 
where Rp is a length scale for disorder and 7' gives the 
strength of the disorder. We took Rp = 0.1 in units of 
the inter-particle separation a, and varied 7/7' in the 
range [0.001:1000]. Without disorder, the ground state 
is a crystal. As the strength of pinning is increased, this 
crystal fragments into smaller domains [15] in a polycrys- 
talline arrangement in which the typical domain size is 
controlled by the pinning strength [16| . We can then cal- 
culate n{B) and its moments by constructing a histogram 



of local field values [17[. In Fig. 4, we compare field con- 



tours obtained for a perfect triangular lattice (a), and 
a disordered arrangement with the attributes discussed 
above (b) . The disordered case shows the absence of long 
range order in both vortex position and field value, while 
maintaining a marked tendency towards local triangular 
coordination. Explicit calculations of a from a distribu- 




FIG. 4: (color online) Contours of constant magnetic field 
strength B(r) computed for (a) a perfect Abrikosov lattice 
at B=0.5 T and for (b) a disordered state with a transla- 
tional correlation length of about 4 inter-vortex spacings gen- 
erated using simulated annealing techniques, with 7/7' = 3 
(see text) . In (c) is shown a plot of the structure factor S{q) , 
obtained analytically as described in the text {q is given in 
units of the inverse mean inter-particle spacing). 



tion of vortex lines with such a structure are a very good 
representation of the experimental data in Fig. 2. 

The same general result may be obtained by analytic 
calculations in certain limits. The triplet structure factor 
is conventionally decomposed as 

5^(qi,q2) = 5(<zi)5(g2)^(l-qi-q2l)[l VC^^'nqi,q2)], 

where C'-^-' is the triplet direct correlation function [l^ 
and p is the average density of the liquid [12] • In our ear- 
lier decomposition of S'^'^^(qi, q2), using the convolution 
approximation, we assumed C'^-* = 0, an approximation 
which is clearly inadequate here [l8|. To progress be- 
yond the convolution approximation, we assume that the 
lines are rigid over the scale of the penetration depth 
and we model translational correlations deep in the vor- 
tex glass phase in terms of correlations in an equilib- 
rium fluid [l3], usin g an accurate analytic approxima- 
tion for C'^'^^(qi, q2) [l^]- We use the pair correlations 
of the jKo{r/X) potential appropriate to rigid flux lines 
to describe local structure in the glass phase. These are 
obtained from solutions to the self-consistent equations 
of classical liquid state theory in the hyper-netted chain 
approximation [13]. Our model S{q) (Fig. 4(c)) is de- 
rived from liquid state computations at a single density 
(A/a ^ 4.0, with a the mean inter-particle spacing, cor- 
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responding to a field of 0.5 T and A ~ 2800 A) using a 
value for the coupling constant T = j/ksT of 50. We 
assume that the local structure of the glass as captured 
in S{q) is not altered substantially as the field is varied, 
once all length scales are expressed in units of the mean 
inter-particle spacing a, an assumption which should be 
valid in the limit where A 3> a. Using Eqs. 1 and 2, 
we calculate a, illustrated in Fig. 2 as the sequence of 
connected blue squares. Note that as the field value is 
increased, the third moment changes sign, with a sat- 
urating at a value of about -0.6, close to the value in 
the experimental data. This relatively simple analytical 
model thus enables a robust description of the effects of 
three-body correlations in the vortex glass phase. 

Similar behavior is also seen in the highly anisotropic 
superconductor Bi2Sr2CaCu208+5 over a range of dop- 
ing [sl, fol , suggesting that system-specific interpretations 
of the negative third moment are unhkely. In ref. 
a transition from a triangular to a square vortex lattice 
was observed in the more highly doped system. We find 
no evidence for this at the field values we probe. At 
high fields our /iSR data do not show the tail in n(B) 
on the high field side which would signal a crystalline 
arrangement, whether square or triangular. Magnetic 
order can coexist with superconductivity in LSCO over 
a restricted doping range. We have characterized the 
sample with neutrons and with longitudinal /x SR, find- 
ing magnetic signatures only below about 4 K [24]. We 
therefore restrict ourselves here to temperatures above 
this value where any magnetic fiuctuations, should they 
exist, lie well outside the muon time window, and cannot 
contribute to the depolarization. 

The sign of the third moment reflects the competi- 
tion between the (positive) contributions from the vortex 
cores, which yield the positive tail of n{B), and (nega- 
tive) contributions from field values at the centres of the 
triangles formed locally by the vortices and associated 
with the minima of n{B). Structures with strong lo- 
cal triplet correlations but no long-range order protect 
both these contributions, but subtly enhance the nega- 
tive ones, since now the positions of particles can fiuctu- 
ate relative to each other, unlike in the perfect crystal, 
while retaining a strong tendency to local triangular co- 
ordination (see Fig. 4(b)) as in the crystal. The sign 
changes arises when the negative contributions to the 
integral (Eq. 1) overwhelm the positive contributions. 
The integrand of Eq. 1 varies strongly as a function of q 
in the vicinity of the first peak of S{q) and below. For 
larger q, form factor cut-offs set in, while at smaller q, the 
integrand becomes negative, due to the generally larg e 
and negative value of C'-^-' in this region (see Ref. [l9|). 
The resultant sign depends on the location of the first 
peak of S (q) , which is itself determined by the magnetic 
field. Thus, non-trivial three-body correlations arise out 
of C(3) in a disordered system (which, unlike a crystal, 
has contributions to the integral at smaller q than the 



first Bragg peak i.e. at length scales which are a little 
larger than nearest neighbor vortex spacings) . Although 
the change in sign of the third moment at high fields con- 
firms the existence of non-trivial three-body correlations, 
it does not indicate the formation of a new vortex state. 
Instead, this observation supports our simple theoretical 
and computational models of structure and correlations 
in the vortex glass state, from which such a change in sign 
follows naturally as the magnetic field is varied, without 
the requirement of a phase transition. 

In conclusion, this paper describes an unusual exper- 
imental consequence of many-particle correlations in a 
magnetic flux line system, showing how three-body corre- 
lations are responsible for negative third moments in the 
field distributions associated with glassy phases of vortex 
lines. Our results motivate and validate the use of sim- 
ple analytic approximations to describe three-body cor- 
relations in bulk disordered systems, an approach which 
should find wider application in areas outside the field of 
superconductivity. 
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